✓s\ 


NASA Technical Memorandum 88295 


* 


Accuracies of Southwell and Force/Stiffness 
Methods in the Prediction of Buckling Strength 
of Hypersonic Aircraft Wing Tubular Panels 


William L. Ko 


(BASA— TB-88295) ACCBBACIES Cl EC01HBE1L AMD M88-1709Q 

FCBCE/S1IFFMKES BE1BCES I* ILF IBEDICTICK OF 
FECK LI KG ST EE MG IB OF HSPEBSCA1C AIBCBAFT 

S2KG 1UEULAE f ABELS (KASA) 28 p CSCL 20K Onclas 

G3/39 0123305 


November 1987 


f 


NASA 

National Aeronautics and 
Space Administration 


NASA Technical Memorandum 88295 


Accuracies of Southwell and Force/Stiffness 
Methods in the Prediction of Buckling Strength 
of Hypersonic Aircraft Wing Tubular Panels 


William L. Ko 

Ames Research Center, Dryden Flight Research Facility, Edwards, California 


l\l/\SA 

National Aeronautics and 
Space Administration 
Ames Research Center 
Dryden Flight Research Facility 
Edwards, California 93523-5000 




SUMMARY 


Accuracies of the Southwell method and the force/stiffness (F/S) method were 
examined when the methods were used in the predictions of buckling loads of hyper- 
sonic aircraft wing tubular panels, based on nondestructive buckling test data. 
Various factors affecting the accuracies of the two methods were discussed. Effects 
of load cutoff point in the nondestructive buckling tests on the accuracies of the 
two methods were discussed in great detail. For the tubular panels under pure com- 
pression, the F/S method was found to give more accurate buckling load predictions 
than the Southwell method, which excessively overpredicts the buckling load. It was 
found that the Southwell method required a higher load cutoff point, as compared with 
the F/S method. In using the F/S method for predicting the buckling load of tubular 
panels under pure compression, the load cutoff point of approximately 50 percent of 
the critical load could give reasonably accurate predictions. 


INTRODUCTION 


The accurate prediction of buckling loads of structural components, based on the 
nondestructive buckling test data, is generally a difficult problem. The graphical 
solution (buckling load prediction) is the test-data plotting schemes, test-data 
extrapolation schemes, and data fitting schemes. The well-known graphical method of 
predicting buckling loads is the Southwell method (or Southwell plot, refs. 1 to 5). 
In the Southwell plot, the compliance (that is, deflection/load) is plotted against 
deflection, and the buckling load is determined from the inverse slope of the plot. 
The Southwell method has been successful in predicting the classical buckling of 
simple structures such as columns and plates (see fig. 1, ref. 6). For complex 
structures exhibiting complex buckling behavior (for example, local instabilities 
and plasticity effect), the Southwell plot may not be a straight line, and therefore 
no discernable slope may be obtained for accurately determining the buckling loads 
( ref . 7 ) . 

The alternate method is the so-called force/stiffness (F/S) method proposed by 
Jones and Greene (ref. 8). In this method, the stiffness (such as force/normalized 
strain) is plotted against the force, and the buckling load is determined from the 
intersection of the limit strain line (the line with slope of unity in the F/S plot) 
and the extrapolation of the curve that fits the data points of the nondestructive 
buckling tests. This method eliminates the concerns about the linearity of the plot. 
The F/S method has been used extensively by Ko, Shideler, and Fields (ref. 9) and by 
Hedges and Greene (ref. 10) in the predictions of the buckling strength of hyper- 
sonic aircraft wing tubular and beaded panels under combined loadings. The F/S 
method seems to give satisfactory buckling load predictions, however, a critical 
review of the accuracy of this method is needed. This report presents the critical 
review of the accuracies of the Southwell method and the F/S method when they are 
applied to the predictions of buckling strength of hypersonic aircraft tubular pan- 
els (structures of complex geometry). 



NOMENCLATURE 


a distance between two adjacent tubes, in 

Bj (j « 1, 2, 3) coefficients of second-order function for least-squares data 

fitting 

C ± (i « 1, 2, • . . , 6) coefficients of third-order function for least-squares 

data fitting 

CSG capacitance strain gage 

D generalized strain variable % 

DT displacement transducer 

E Young's modulus, lbf/in 2 

F applied load, lbf 

F* maximum applied load in nondestructive buckling tests, lbf 

F/S force/stiffness 

f width of panel flat region, in 

k extrapolation factor, F cr /F* 

Z length of tubular panel, in 

m exponent in the expression of D 

N x panel compression load, lbf/in 

N X y panel shear load, lbf/in 

p lateral pressure, lbf/in 2 

R radius of tubular wall, in 

RSG rosette strain gage 

SG axial strain gage 

s length of circular arc element of panel tube cross section * 

T temperature, °F 

t thickness of tubular wall, in 
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t effective thickness of tubular panel, in 

w width of tubular panel, in 

a half angle of circular bead arc 

strain due to bending, in/in 
e c strain due to axial compression, in/in 

y strain due to shear, in/in 

6 lateral displacement at center of tubular panel, in 

[ ] cr critical value of [ ] 


TUBULAR PANEL 


Ko, Shideler, and Fields conducted extensive nondestructive buckling tests 
(ref* 9) of tubular panels made of two formed Rene 41 alloy sheets seam-welded 
together to form five flat regions (double sheets) and four noncircular tubular 
regions (fig. 1). At the end of these tests, panels 1 and 3 were loaded up to buck- 
ling failure under pure compression. Because the actual buckling loads of the tubu- 
lar panel under pure compression are known, it is possible to examine the accuracy of 
the Southwell method and F/S method in predicting the buckling strength of the tubu- 
lar panels under pure compression. The test data of panel 3 will be used in this 
accuracy study. Figure 2 shows the locations of the strain gages and the 11 lateral 
displacement transducers (DT) instrumented on panel 3 for pure compression tests, 
reported in reference 9. Figure 3 shows the buckled sites at the rear surfaces of 
tubes 1, 2, and 4 of test panel 3. 

The data obtained from rosette strain gages RSG622, RSG623, and RSG633? axial 
strain gages SG516 and SG517; and the displacement transducer DT6 will be used in 
the study of accuracies of the af orementioned two methods of predicting buckling 
loads. That instrumentation was located in the center region of the panel and was 
relatively close to the buckled sites on tubes 1 and 2. Figure 4 shows the out-of- 
plane displacement of panel 3 at the time of buckling failure under pure compression 
(failure load F cr - 41,051 lb, ref. 9). Notice that the deformation of the panel is 

a fundamental mode deformation with slight distortion near the buckled zone. 


SOUTHWELL PLOT 


The Southwell method has been successful in the prediction of classical bucking 
of simple structures such as columns and plates (refs. 1 to 5). For complex struc- 
tures, such as a tubular panel with complex buckling behavior, the dependability 
of the Southwell method must be carefully examined. Using the results of the 
destructive buckling tests reported in reference 9, a Southwell plot was constructed 
for a tubular panel (panel 3), and is shown in figure 5. In the figure, F is the 
compressive load, and 6 is the lateral displacement at the center of the panel. 
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It is seen that the data points lying in the region 0.43 < F/F cr < 0.95 (F cr being 

the critical load at buckling) form a nice straight line. However, the inverse slope 
of this straight line yields the critical load of F cr = 55,184 lb which is 34 percent 

higher than the actual buckling load of F cr = 41,051 lb. Also, one notices that the 
data points lying within the region 0.96 < F/F cr < 1.00 form another straight line 
with a slightly steeper slope, yielding the critical load of F cr = 53,846 lb. This 
value is approximately 31 percent over the prediction of the actual buckling load. 

The second Southwell plot for tubular panel 3 shown in figure 6 is the plot 
of 6/N x as a function of 6, where N x is the effective panel compressive load calcu- 
lated from the outputs of the two axial strain gages SG516 and SG517 located at the 
outer surfaces of tube 2 (see figs. 3, 7), namely: 

N x = | SG516 + SG517 | (1) 


where E is the Young modulus and t ( =0.037 in) is the effective thickness of the 
tubular panel. Similar to figure 5, the data points in figure 6 also form a nice 
straight line in the region 0.43 < F/F cr < 0.95. The inverse slope of this straight 

line gives the critical panel load of (N x ) cr = 2687 lb/in, which is 26 percent higher 
than the actual panel buckling load of (N x ) cr = 2138 lb/in ( * F cr /w = 41 ,051/19.2) , 

where w is the panel width (see fig. 1). Also, the data points lying in the region 
0.96 < F/F cr < 1.00 form another straight line with a steeper slope, giving (N x ) cr = 

2553 lb/in, which is 19 percent above the measured value. 

The final Southwell plot shown in figure 7 is the plot of 6/N x as a function of 
6, where the effective panel load N x is calculated from the output of the rosette 

strain gage leg RSG622 located at the flat region of the center of the panel (see 
fig. 2), namely: 


N x = Et | RSG622 | (2) 

The plot, with the exception of a small strain region, is practically bilinear in 
shape. The inverse slope of the straight line fitting the data points lying in the 
region 0.43 < F/F cr < 0.95 gives the critical panel load of (N x ) cr = 2769 lb/in, 

which is a 30 percent overprediction. The steeper straight line fitting the data 
points in the region 0.96 < F/F cr < 1.00 yield (N x ) cr = 2400 lb/in gave a buckling 

load approximately 12 percent higher than the actual panel buckling load of (N x ) cr = 
2138 lb/in. 

From the preceeding three Southwell plots for the tubular panel (see figs. 5 to 
7), it is noticed that the plots are practically bilinear in shape (excluding the 
small strain region below F/F cr = 0.43), and that the F/F C r = 0.95 point lies almost 
at the middle of the plot. This means that the high loading data points in the 
region 0.95 < F/F cr < 1.00 spread over one-half span of the plot. Thus, in order 

to obtain the straight line zone for most accurately determining the slope, the non- 
destructive test has to be carried out up to 90 - 95 percent of the buckling load. 

For the case of the tubular panel, the Southwell method excessively overpredicts 


% 
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the buckling load (either using the slope for the region 0.43 < F/F cr < 0.95, or for 
the region 0.96 < F/F cr < 1.00). 

Figure 8 shows the Southwell plot for a square plate containing a central cir- 
cular hole subjected to pure compression (ref. 7). For this plot, the ratio of hole 
diameter to plate length is 0.1. It is clear that the Southwell plot for this par- 
ticular perforated plate is nonlinear, having no discernable slope for determining 
the buckling load. 

During the nondestructive buckling tests, the Southwell plot may be displayed, 
say, on the screen of a cathode ray terminal (CRT) for observation of the growth of 
the plot. The test can be terminated if enough data points have formed a straight 
line. However, it is impossible to know the loading ratio F/F cr of the test cutoff 
point because the location of the possible buckling point is unknown. As will be 
seen shortly, this problem does not occur in the F/S method. 

FORCE/STIFFNESS METHOD 

The force/stiffness (F/S) method, which is illustrated in figure 9, was devel- 
oped by Jones and Greene (ref. 8) primarily for predicting the local buckling failure 
of complex structures (for example, beaded and tubular panels for hypersonic aircraft 
wings) under combined loading. The method was used extensively by Ko, Shideler, and 
Fields (ref. 9), and Hedges and Greene (ref. 10) in the nondestructive buckling 
strength predictions of beaded and tubular panels for a hypersonic aircraft wing. 

In this method, force F (which can be any particular load of applied load set) is 
plotted against stiffness F/D, where D is the generalized strain variable (refs. 9 
to 10), and the buckling load is obtained at the intersection of the extrapolation of 
the curve fitting the nondestructive buckling test data points and the limit strain 
line (the line with slope of unity for which D - D cr - 1, see fig. 9). The general- 
ized strain variable D is defined as (ref. 8) 



where e c , e^, and y are respectively the strains in compression, bending, and shear, 

[ )cr denotes critical values, and the exponent m is an empirically determined 
constant. The values for e c , and y are obtained from the outputs of strain gage 

measurements, and (e c )cr/ ( e b)cr* and Ycr must be calculated analytically using the 

buckling equations pertaining to panel local geometry (for instance, circular cylin- 
drical panel for the tubular panel, refs. 8 to 10). At the critical local buckling 
strain state, equation (3) is set to unity, namely: 

D = D cr = 1 (4) 

which represents a critical strain interaction surface, as shown in figure 10. 

Figure 1 1 shows the F/S plot of the same set of pure compression-destructive 
buckling data used in the earlier Southwell plots (see figs. 5 to 7). Unlike the 
Southwell plots, the data point of F/F cr 38 0.95 falls very close to the limit strain 
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line (D cr = 1) since the abscissa is F instead of 6. This means that when the F/S 

method is used, the nondestructive buckling test can be cut off at a relatively lower 
loading level compared with the Southwell method (see fig. 5). The fact that the 
actual buckling point is very near the limit strain line implies the accuracy of 
the F/S plot. Notice that the data points in the region 0.3 < F/F cr < 0.95 form a 

smooth curve , and the curve that fits those data points intersects with the limit 
strain line at a point yielding a predicted buckling load almost equal to the actual 
buckling load. 

Figure 12 shows the F/S plot (taken from ref. 9) for combined compression and 
shear loading N xy /N x - 1.22. When there is no lateral pressure acting, the F/S plot, 

except for the small load region, is usually strongly convex downward, allowing the 
least-squares data fitting to give relatively accurate predictions of the buckling ^ 

load. When there is lateral pressure in addition to the combined compression and 
shear loading, the F/S plot usually turns out to be convex upward. This is shown 

in figure 13 (taken from ref. 9), for which N xy /N x = 3.21 and p * 0.75 lb/in 2 * For 
a plot shape, like the one shown in figure 13, sometimes more accurate buckling load 
prediction could be obtained by curve-fitting only the higher load region and ignor- 
ing the lower load region, which is less important. 


Accuracy of Force/Stiffness Method 

The accuracy of the force/stiffness (F/S) method depends mainly on the follow- 
ing factors: 

Test cutoff point . — The higher the load at which the nondestructive buckling 
test is cut off, the better the prediction of buckling load because the range of 
extrapolation of the data-fitting curve is shorter (or lower extrapolation factor k « 
F cr /F*, see fig. 9). 

Mathematical function used in curve fitting . — The predicted buckling load will 
vary with the mathematical function chosen for the least-squares fitting of the data 
points. For combined compression, shear and lateral pressure loading, the F/S plot 
is always convex upward (see fig. 13). The least-squares fitting curve, based on a 
certain mathematical function, may tend to bend upward outside the data point range 
and may cause overprediction of the buckling load, or may not intersect with the 
limit strain line. 

Range of curve fitting . — The data points in the lower load region are less 
important than the data points in the higher load region. The accuracy of the 
buckling load prediction could sometimes be improved by curve-fitting only the 
data points in the higher load region and ignoring the data points in the lower 
load region. 

Exponent m . — The value of the exponent m in the expression of D will affect 
the value of D. The value of m for most types of panels was found to be nearly 2 
(ref. 10). Thus, the effect of m on D is minimal. 

Accuracy of (e c ) cr , Ub^cr* Ter* “ The value of D is affected by the accuracy of 
the calculated (e c )cr# ( e b)cr# an<i Ycr* Since the buckling theories used in the 
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calculations of (e c ) cr , (e b ) cr , and y C r are well established, the effect of the values 
of (e c ) cr * ( e b^cr# and Ycr on the vaiu e of D or F/D could be very small. 

Density of strain gage sites . — In using the F/S method, the test structure is 
instrumented with strain gages at different sites for monitoring local buckling. The 
greater the number of strain gage measurement sites, the better the chance to have 
some strain gage sites located near the potential buckle sites, and, therefore, a 
greater downward-bend region of the F/S plot may be obtained for more accurate pre- 
diction of the buckling load. Namely, the wider the downward-bend region is, the 
more accurately the data-fitting curve can give its extrapolation for determining the 
buckling load. 
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Test Cutoff Point 

As previously mentioned, the most critical factor affecting the accuracy of the 
F/S plot is the location of the test cutoff point (for example, how far the test 
cutoff point is from the limit-strain line). The data presented in figure 11 may be 
used for the study of the effect of the location (or load level) of the test cutoff 
point on the accuracy of the F/S method. The load cutoff points selected for the 
preceding study were F/F cr = 0.5, 0.6, 0.7, 0.8, 0.9, 0.95, and the following two 

types of mathematical functions were used in the least-squares data fittings: 

F 

1. Second-order function: — = Bi + B 2 F + B 3 F 2 

2. Third-order function: 

F 1 + CjF + C2F 2 

D c 3 + C4F + C5F 2 + C 6 F 3 

where Bj (j * 1, 2, 3) and Ci (i = 1, 2, ... 6 ) are constants determined from the 
least-squares fittings of the test data. 

Figures 14 and 15 respectively show the F/S plots of test data of figure 11 
fitted by the second-order and third-order functions up to different load cutoff 
points. In the least-squares data fittings, the insignificant data points in the 
region F/F cr < 0.3 were neglected. The least-squares fitting for F/F cr = 0.8 using 

the third-order function (equation ( 6 )) gave poor extrapolation because the discon- 
tinuity (region of sudden slope change) fell near the force cutoff point F/F cr = 0.8. 

For other values of F/F C r* the discontinuity did not occur within the test data 
region 0.3 < F/F cr < 1.0, and therefore, it did not affect the buckling load predic- 
tions. The third-order function (equation ( 6 )) was originally constructed for 
fitting the nondestructive buckling test data exhibiting sudden slope change (see 
figs. 9 and 12). For fitting the test data having a smooth slope change, equa- 
tion ( 6 ) may not be ideal for certain values of F/F cr . Overall it appears that the 

second-order function will give the most consistent prediction of buckling strength 
of tubular panels under pure compression. Table 1 summarizes different values of 
pure-compression buckling loads predicted from the plots in figures 14 and 15. The 
table also lists the error involved in each buckling load prediction. For the load 
cutoff points F/F cr = 0.5, 0.9, 0.95, the third-order function gave more accurate 


(5) 


(6) 
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predictions. However in the region 0.6 < F/F cr < 0.8, the reverse is true. The data 

given in Table 1 suggest that when the F/S method is used in the prediction of buck- 
ling of tubular panels under pure compression, the nondestructive test may be termin- 
ated at a relatively low load cutoff point, say F/F cr * 0.5. 


CONCLUSIONS 


Accuracies of the Southwell method and the F/S method were critically reviewed 
for the case when they were applied to the prediction of buckling loads of hypersonic 
aircraft wing tubular panels, based on the nondestructive buckling test data. Various 
factors that affect the accuracies of the two methods were discussed in great detail. 
For the case of tubular panels, the F/S method gave more accurate buckling load pre- 
diction than the Southwell method, which excessively overpredicts the buckling load. 
For pure compression, the F/S method could give reasonably accurate buckling load 
predictions, even with the load cutoff point as low as F/F cr =* 0.5. The Southwell 
method requires a much higher load cutoff point. 


r 
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TABLE 1. ACCURACIES OF BUCKLING LOADS IN PURE COMPRESSION PREDICTED 
FROM FORCE/STIFFNESS METHOD USING DIFFERENT LOAD CUTOFF POINTS AND 
LEAST-SQUARES DATA FITTING FUNCTIONS. 


Load 

cutoff point 
F/F cr 

Predicted buckling load, 
F cr , lb 

Error, percent 

Second-order 

function 

Third-order 

function 

Second-order 

function 

Third-order 

function 

0.5 

43,350 

39,770 

5.60 

3.12 

0.6 

44,660 

45,840 

8.79 

11.67 

0.7 

42,930 

45,430 

4.58 

10.67 

0.8 

42,690 

46, 1 10 a 

3.99 

12.32 a 

0.9 

42, 160 

40,580 

2.70 

1.15 

0.95 

41,950 

41,170 

2.19 

0.29 

1.0 

41,051 b 

— — 

0 

0 


a Poor extrapolation curve. 
^Actual buckling load. 


/ 


- _ 2t (» + f) 
a 


= 0.037 


7172 


Dimensions in inches . 
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Figure 3. Buckled tubular panel 3 , rear side s after 
pure compression test . 
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Figure 7 . Southwell plot for tubular panel 3 under pure compree8ion. N x calcu- 
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Figure 9. Force/stiffneae plot for local buckling . 
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Figure 14. Least-squares fittings of different regions of pure compression buckling test data of 
tubular panel using second-order function. 
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